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Abstract. We further investigate the runaway instability 
of relativistic tori around black holes assuming a power law 
distribution of the specific angular momentum. We neglect 
the self-gravity of the torus; thus the gravitational field of 
the system is determined by the central rotating black hole 
alone (the Kerr geometry). We find that for any positive 
(i.e. consistent with the local Rayleigh stability condition) 
power law index of the angular momentum distribution, 
the tori are always runaway stable for a sufficiently large 
spin of the black hole. However, before concluding that 
some (astrophysically realistic) tori could indeed be run- 
away stable, one should include, in full general relativity, 
the destabilizing effect of self-gravity. 
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1. Introduction 

Toroidal fluid conflgurations around black holes, known 
as thick accretion disks, have been suggested as models of 
quasars, other active galactic nuclei, and some X-ray bina- 
ries. Recently they have also assumed an important role in 
cosmological models of 7-ray bursts based on the merging 
of two neutron stars. A rather simple mathematical class 
of such configurations has been introduced by Fishbone & 
Moncrief (1976), Fishbone (1977), and fully described an- 
alytically by the Warsaw group: Abramowicz, Jaroszynski 
& Sikora (1978), Kozlowski, Jaroszynski & Abramowicz 
(1978), Jaroszynski, Abramowicz & Paczynski (1980), and 
Paczynski & Wiita (1980). In the last paper, a very prac- 
tical and accurate Newtonian model for the gravitational 
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field of a non-rotating black hole, known as the Paczyhski- 
Wiita potential, was introduced. 

The question of stability of such configurations at- 
tracted considerable attention, because it was recognized 
that some types of instability could have very direct, and 
quite interesting, astrophysical consequences. Obviously, 
for the same reason as in the Newton theory, the essen- 
tially local Rayleigh criterion for dynamical stability with 
respect to axially symmetric local perturbations demands 
that the specific angular momentum / should increase with 
the distance R from the axis of rotation (Seguin 1975). 
Even the Rayleigh-stable tori with dl/dR > are dy- 
namically unstable in the presence of a weak magnetic 
field (Balbus & Hawley 1991). The Balbus-Hawley insta- 
bility does not destroy the large scale structure of tori, but 
instead drives a local turbulence which induces viscosity 
that is needed for accretion to occur. Much more threaten- 
ing for the global structure was the important discovery by 
Papaloizou & Pringle (1987) that all non-accreting tori are 
unstable to global non-axisymmetric perturbations. The 
consequences of this brilliant work have been studied by 
numerous authors, and the present view is that even a 
very modest mass loss due to accretion may be sufficient 
to stabilize the Papaloizou-Pringle modes (Blaes 1988). 
Also, it has been shown that the self-gravity of the disk 
has a stabilizing effect (Goodman & Narayan 1988). 



2. The runaway instability 

A general feature of fluid distribution around a black hole, 
is the existence of a self-crossing equipotential surface, 
similar to the Roche lobe in binary systems. The circle 
at which it crosses itself is called the cusp. The material 
inside the Roche lobe is gravitationally bound and only in 
such a case is equilibrium possible. Those configurations 
that just fill their Roche lobes are marginally bound and 
thus could be called critical. The location of the cusp and 
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the shape of the relativistic Roche lobe is determined by 

the combined gravitational potential of the black hole and 
the torus, and by the centrifugal potential due to the disc 
rotation. 

For critical tori, when a small amount of material is ac- 
creting into the hole, a natural question arises (Abramow- 
icz, Calvani & Nobili 1983): is the torus stable with respect 
to the transfer of mass through its cusp- like inner edge? To 
answer this question, one must determine the new position 
of the cusp, the shape of the Roche lobe, and the new equi- 
librium for the torus. If the relativistic Roche lobe shrinks 
sufficiently enough, matter which was bounded before the 
mass transfer will become unbounded, falling catastroph- 
ically into the hole on a dynamical time scale. This is the 
runaway instability. A full answer to the question about 
the conditions for the runaway instability to occur is far 
more difficult than was originally imagined, and is still 
unknown today. The claim by Abramowicz et al. (1983) 
that sufficiently massive I = const tori are all runaway un- 
stable, was based on an approximate model in which the 
gravitational field of the central black hole was modeled 
by the Paczyhski-Wiita potential. The self-gravity of the 
torus was included using Newtonian gravity. However, us- 
ing the Kerr metric for the black hole, but ignoring the 
self-gravity of the torus, Wilson (1984) demonstrated that 
the non self- gravitating tori do not suffer from the run- 
away instability. The reason is that accreting matter de- 
livers not only mass but also angular momentum to the 
black hole in such a way that the ratio a/M (dimension- 
less specific angular momentum of the black hole) remains 
approximatively constant thereby leaving unchanged the 
structure of the space time. By applying another approach 
(a Newtonian model of a rotating black hole-ring sys- 
tem) Khanna & Chakrabarti (1992) confirmed however 
that self-gravitating rings are runaway unstable. 

It is obvious that the situation is complex because the 
self-gravity makes tori runaway unstable, while rotation of 
the central black hole stabilizes them. Thus, both effects 
are important and neither can be neglected. To clarify the 
issue, Nishida et al. (1997), using numerically constructed 
full general-relativistic models of self-gravitating tori with 
I = const found that those that are sufficiently massive 
(the ratio = M^/MBH > 0.1 of the disc mass 
to the central black-hole mass M) are runaway unstable 
always. These authors pointed out that the runaway insta- 
bility introduces an acute difficulty for the recently consid- 
ered models for 7-ray bursts based on the merging of two 
neutron stars by making the lifetime of the torus (which 
is a product of the merging) too short to provide enough 
energy to power the burst. 

Most recently, however, Daigne & Mochkovitch (1997) 
made an important discovery that a non-zero gradient of 
the angular momentum distribution inside the torus has 
a strong stabilizing effect. They used models of non self- 
gravitating tori and the Paczynski-Wiita potential to de- 
scribe the black hole. In this case, after the mass transfer. 



the angular momentum content in the innermost part of 

the resulting configuration is higher and therefore provides 
a stronger gravitational barrier preventing catastrophic 
runaway accretion. 

To verify that such a stabilizing effect is also present 
in full general relativistic regime, one should study the 
stability of equilibrium configurations constructed self- 
consistently, in a similar way as in Lanza (1992) and 
Nishida & Eriguchi (1994) but for a non-constant dis- 
tribution of angular momentum. Constructing general- 
relativistic equilibrium configurations with an arbitrary 
Eulerian rotation law is not particularly difficult. Ex- 
plicit analytic models with general Eulerian radial dis- 
tribution of angular momentum, I = 1{R) (for exam- 
ple with a power-law I (x R"^) have been constructed by 
Jaroszyhski et al. (1978), and later used by Kuwahara 
(1988), Chakrabarti (1991), and others. If the equation of 
state is fixed, each particular model in this class is de- 
termined by the ratio of the disc mass to the 
central black-hole mass M and by the location of its in- 
ner radius. However, since we want to study the stability 
of a torus that is accreting axisymmetrically, it is neces- 
sary to mimic the process by means of a quasi-stationary 
sequence of equilibria along which the Lagrangian distri- 
bution of angular momentum per baryon j{m) is con- 
served. From a numerical point of view this is a diffi- 
cult task. Already in Newtonian theory finding equilib- 
ria with a given Lagrangian angular momentum distribu- 
tion demands a non trivial effort (Ostriker & Marck 1968, 
Daigne & Mochkovitch 1997). The self-consistent iteration 
between the Einstein field equations and the equation de- 
scribing the equilibrium of matter is complicated by the 
fact that the assumed distribution of the angular momen- 
tum depends on the Lagrangian coordinate m, whereas the 
field equations are expressed in terms of Eulerian coordi- 
nates (e.g. Kozlowski et al. 1978). (Of course, this prob- 
lem would disappear if we had spherical symmetry, since 
in this case one could easily rewrite Einstein equations in 
terms of m.) Moreover, in general relativity, the surfaces 
where the specific angular momentum is constant (von 
Zeipel surfaces) do not have a trivial cylindrical structure 
R = const, as in Newtonian theory, but must be deter- 
mined iteratively. Their shape is more complicated, given 
as an implicit function TZ{R, 0) = const, and it is changed 
during accretion. This fact is essential for construction of 
the tori. It is convenient therefore to explore such compli- 
cated iterative techniques first in a simpler approximation 
before solving the full problem. 

In this Note we have solved a particular technical prob- 
lem that is necessary to obtain the full solution, namely 
we established a method of finding non self-gravitating 
equilibria with a given Lagrangian distribution of angular 
momentum in general relativity. Apart from a straightfor- 
ward replacement of Newtonian formulae by correspond- 
ing relativistic expressions, the major difference and diffi- 
culty in studying relativistic test-fluid tori which we have 
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faced in this work follows from the fact that the structure 
of von Zeipel surfaces is not given a-priori and must be 
determined iteratively. 

3. Numerical procedure 

We assume that the spacetime is axially symmetric, sta- 
tionary, and described by the Kerr metric g^i, in Boycr- 
Lindquist coordinates {t, r, 0, 0} (geometric units with 
c — G = 1 in usual notation). The metric is parame- 
terized by the black-hole mass M and its specific angular 
momentum a/M. The torus is described by a relativistic 
polytrope which relates the pressure P to the density of 
material p: P = np^ with k and 7 constant (Tooper 1965). 
The corresponding energy density e = p -t- P/(7 — 1) is re- 
lated to the energy per unit mass £ — {P + e)p~^ut and 
to Lagrangian angular-momentum density j = £1 {ut de- 
notes the time-component of the four- velocity of the fluid; 
cf. Kozlowski et al. 1978). 
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Fig. 1. Illustration of the iterative procedure for constructing a 
model of the torus with a prescribed mass. See text for details. 



First, an initial critical configuration is constructed, 
given the parameters M^, M, a, k, and 7. The starting 
angular-momentum distribution is specified as a power- 
law I oc TZ'^ with q > constant. The model is then de- 
termined by the above-given parameters, in particular the 
inner and the outer edges, and the centre of the torus 
where the pressure is maximum. The matter distribution 
inside the disc is calculated and expressed in terms of a 
spline-interpolated function of the von Zeipel coordinate, 
m{TV), where m(7?.(rin)) = and normalisation guarantees 
m(7^(rout)) = 1 at the outer edge. Also the Lagrangian 
angular-momentum density j is expressed as a normalised 
spline- interpolated function: 2 = 1(7^), < Z < 1. This 
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Fig. 2. Shading indicates the region in space of parameters 
{q,a/M) in which the tori filling their Roche lobes are run- 
away unstable with respect to small overflow. Dark-shaded area 
corresponds to the initial mass of the torus Af^ = 0.36Mq, 
light-shaded area refers to AI^ — O.I8M0. Black-hole rota- 
tion has apparently a stabilizing effect. The horizontal line 
a.t q — 0.07 indicates where the stability boundary is located 
in the corresponding pseudo-Newtonian analysis of Daigne & 
Mochkovitch (1997). See text for details. 



initial configuration is computed with a pre-determined 
accuracy of . Figure |^ illustrates the above-described 
iterative procedure. Each bundle of the curves corresponds 
to a set of surfaces with the inner edge fixed and increasing 
resolution (the effect of improving resolution is clearly vis- 
ible). Construction of the models is terminated when the 
mass of the torus is determined with accuracy better than 
a pre-determined value, typicaly 0.1%. Location of the in- 
ner edge is then changed (another bundle of the curves) 
and the code iterates again, producing another configura- 
tion with a different mass of the torus. The whole proce- 
dure is repeated (five times in this example) with the in- 
ner edge bounded between the marginally bound and the 
marginally stable orbits. Iterations are terminated when 
the torus has a right mass. The final configuration is in- 
dicated by shading. 

Now we assume that a small amount of material con- 
tained near the cusp, between rin and rin + Sr-m, is trans- 
ferred throughout the cusp to the black hole. The mass 
ratio rrij^ = M-^/M is decreased by 5m-^, typicaly of the 
order of one percent, while the angular-momentum pa- 
rameter of the black hole is changed by a corresponding 
value, 5m^l{ri^). The distribution of j is conserved dur- 
ing the process of accretion: j'{m') = j(m) (primes denote 
new quantities, after the mass transfer). The next iteration 
starts with a guess for the new function Iq ~ j'/£- The 
inner edge of the new critical configuration is iterated un- 
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til a configuration with the required rn^, is found; at each 
iteration I'q is also updated, so that the final r(,-, and I' (TV) 
are known only when convergence has been reached. This 
updating induces a corresponding change of the von Zeipel 
surfaces, TZ(r, 9) — > TZ'{r, 9) which provides however only a 
minor modification and is visible near the inner edge. At 
convergence then one can verify whether the inner edge 
of the new critical surface has moved. If + Sri^ > r[^, 
then the new configuration is unstable since the crossing 
equipotential has shrunk. 

4. Results and Conclusions 

We have adopted the same values of the starting param- 
eters as Daigne & Mochkovitch (1997): M = 2.44 Mq, 
7 = 4/3, K = 0.5 X 10^^ (cgs), corresponding to the mod- 
els of presumed massive tori which are formed in neu- 
tron star mergers. The mass of the torus was chosen to 
be = 0.36 Mq and M„ = 0.18 M©, Sm„ = 1%. We 
assumed that the torus corotates with the black hole. Fig- 
ure H illustrates the results of the calculation, namely the 
boundary between stable and unstable configurations in 
the plane of parameters (g, a/M) of the initial configura- 
tion. It can be seen that increasing q and increasing a/M 
both stabilize the configuration. This is the main result 
of our present note. The boundary between stability and 
instability is positioned at roughly the same value of q 
as found by Daigne & Mochkovitch (1997) in their work. 
The difference is factor of 2 which cannot be considered 
too large given the approximate nature of the Paczyhski- 
Wiita potential and the difference in the iterative pro- 
cedures. Naturally, Daigne & Mochkovitch could not ex- 
plore the dependence on a/M which does not appear in 
the Paczyhski-Wiita potential they adopted. The horizon- 
tal line in Fig. || indicates where the boundary is located 
in the corresponding analysis of Daigne & Mochkovitch 
with = 0.36 Mq when pseudo-Newtonian radial coor- 
dinate is formally identified with Boyer-Lindquist r in our 
analysis. 

It remains to be discussed how self-gravity of the disc 
affects this criterion of stability of the accretion process, 
and what is the interplay between effects of a rotating 
black hole and a self-gravitating torus within full general 
relativity (work in progress). 
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